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1—5 | Abstract 

. Fractional generalization of an exterior derivative for calculus of variations is defined. 

The Hamilton and Lagrange approaches are considered. Fractional Hamilton and Euler- 
Lagrange equations are derived. Fractional equations of motion are obtained by fractional 
qq ■ variation of Lagrangian and Hamiltonian that have only integer derivatives. 

-xj- : 

PACS numbers: 45.20.-d; 45.20. Jj; 45.10.Hj 

o ■ 

O ■ 1 Introduction 

The theory of derivatives of noninteger order P3 [2] goes back to Leibniz, Liouville, Riemann, 
Grunwald and Letnikov 2]. Derivatives and integrals of fractional order have found many 
applications in recent studies in mechanics and physics. In a fairly short period of time the list 
of such applications becomes long. For example, it includes chaotic dynamics [3111!, mechanics of 
fractal media 00 13 El, quantum mechanics 0QIIJ, physical kinetics [31 EEH EE2 ECU El E3 EE] , 
plasma physics [T71 HHj, astrophysics [T9*j . long-range dissipation [201 [21], mechanics of non- 
Hamiltonian systems [22 [23] , theory of long-range interaction [2H [23 EH] , anomalous diffusion 
and transport theory [31 [23 12H] and many others physical topics. 

In mathematics and theoretical physics, the variational (functional) derivative is a general- 
ization of the usual derivative that arises in the calculus of variations. In a variation instead of 
differentiating a function with respect to a variable, one differentiates a functional with respect 
to a function. In this paper, we consider the fractional generalization of variational (functional) 
exterior derivatives. 

The main results are derived in sections 4.2, 5.2, 6.2, and 6.3. In sections 2, 3, 4.1, 5.1, 6.1, 
brief reviews of fractional derivatives, differential forms, Hamiltonian systems are considered 
to fix notation and provide convenient references. In section 2, a brief review of differential 
forms is considered. In section 3, we consider Hamiltonian and fractional Hamiltonian systems 
|23j . In section 4, we define the fractional variations in Hamilton's approach to describe the 
motion. The fractional generalization of stationary action principle is suggested. In section 
5, we discuss the fractional variations in Lagrange's approach to describe the motion, and the 
fractional generalization of stationary action principle is suggested. In section 6, we consider 
the generalization of action principle to non-Hamiltonian systems. The fractional equations of 
motion with friction are discussed. Finally, a short conclusion is given in section 7. 



1 



2 Fractional Derivatives and Differential Forms 
2.1 Differential Forms 

In this subsection, a brief review of differential forms fl§\ I3*U] is considered to fix notation and 
provide a convenient reference. 



Definition 1. A differential 1-form 

co = F l (x)d,Xi 

is called an exact 1-form in R n if the vector field F l (x) can be presented as 



dV 



where V = V(x) is a continuously differentiate function. 



In this case, the differential form (0) is an exact form u = —dV, where V = V(x) is 
a continuously differentiate function (0-form). Here d is the exterior derivative [21]. The 
exterior derivative of the function V is the 1-form dV = dxi dV/dxi written in a coordinate 
chart (xi, ...,x n ). For the k-form u k and the /-form ui, the exterior derivative obeys the relation 

d(u k A uji) = (du k ) Auji + (-l) k u k A dwj. (3) 

Here k and I are integers. Note that dduj = for any form u. If du = 0, then uj is called a 
closed form. 

In mathematics the concepts of closed and exact forms are defined by the equation 
du = for a given u to be a closed form, and u = dh for an exact form. It is known that 
to be exact is a sufficient condition to be closed. In abstract terms the question of whether 
this is also a necessary condition is a way of detecting topological information by differential 
conditions. 



Proposition 1. If a smooth vector field F = e{F l (x) satisfies the relations 



dxj dxi ^ 
on a contractible open subset X of R n , then (QJ) is the exact form such that 

9V{x) , 

u = dx^ (5) 

dx, t 

Proof. Let us consider the forms (0). The formula for the exterior derivative of (0) is 

1 (3F i 8Fi\ , 

duJ = 2\-dx- ~ ^) jA * 
where A is the wedge product Therefore the condition for u to be closed is (jlj). If 

F % = —dV/dxi, then the implication from 'exact' to 'closed' is a consequence of the per- 
mutability of the second derivatives. For the smooth function V = V(x), the second derivative 
commute, and equation (j3J) holds. 
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2.2 Fractional Differential Forms 

If the partial derivatives in the definition of the exterior derivative 

dxi 

are allowed to assume fractional order, then a fractional exterior derivative is defined [31] by 

d a = {dxiYT>« x , (6) 

Here we use 

D 7 (,) = 1 r d JL (7) 

xJK ' T{m-a)Jo {x-y) a ~ m+1 dy m ' 1 ' 

where a > 0, and m is the first whole number greater than or equal to a. Equation (JZJ) defines 
the Caputo fractional derivatives jHSl 1331 El ESJ of order a > 0. 

Definition 2. A differential 1-form 

u a = F l (x)(dxi) a (8) 
is called an exact fractional form if the the vector field F t (x) can be represented as 

F\x) = -D« V, (9) 
where V = V(x) is a continuously differentiable function, and is a derivative of order a. 

Using (jHJl the exact fractional form can be represented as 

Ua = -d a v = -{dxi) a TiZy. 

Therefore, we have Q. 

Note that equation (jHJ) is a fractional generalization of the differential form ((TJ. Obviously 
that fractional 1-form u a can be closed when the differential 1-form uj = uj\ is not closed. The 
fractional analogue of proposition 1 has the form. 

Proposition 2. If a smooth vector field F = eiF l (x) on a contractible open subset X of R n 
satisfies the relations 

B^ j F i - B^F j = 0, (10) 
then the form is an exact fractional 1-form such that 

w a = -D«V(x), (11) 

where V(x) is a continuous differentiable function and D".V = — F l . 

Proof. This proposition is a corollary of the fractional generalization of the Poincare lemma 
The Poincare lemma is shown j^U E2] to be true for exterior fractional derivative. 
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Note that we can generalize the definition of fractional exterior derivative by the equation 

n 

d a = Y,(dx t rD: : , (12) 

i=l 

where a = (a±,a2, ...,a n ), and consider the fractional differential 1-forms: 

n 

Lu a = J2ui(x){dxi) a '. (13) 
i=i 

In this case, we can derive equations with derivatives of different orders Oj. For simplicity, we 
suppose that all a, = a. 

3 Hamiltonian Systems 

In this section, a brief review of Hamiltonian systems [21] and fractional Hamiltonian systems 
[23] is considered to fix notations and provide a convenient reference. 

3.1 Definition and Properties of Hamiltonian Systems 

Let us use the canonical coordinates (qi, ...,q n ,pi, ...,p n ) in the phase space R 2n . We consider 
a dynamical system that is defined by the equations 

The definition of Hamiltonian systems can be realized in the following form [SHI I2HJ EE| ■ 

Definition 3. The dynamical system fli^| ) on the phase space R 2n , is called a Hamiltonian 
system if 

(5 = GUpi - F*dqi, (15) 

is a closed form, d/3 — 0. A dynamical system is called a non- Hamiltonian system if J73)) is 
non-closed, d/3 0. 

The exterior derivative for the phase space is defined as 

d d 

d = dqi— + dpi — . (16) 
dqi dpi 

Here and later we mean the sum on the repeated index i from 1 to n. 

Proposition 3. // the right-hand sides of equations satisfy the conditions 

dpj dpi ' dqi dpj ' <% % 
then the dynamical system is a Hamiltonian system. 
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Proof. In the canonical coordinates (q,p), the vector fields that define the system have the 
components (G l ,F l ), which are used in equation (|14J) . Let us consider the 1-form (|15jl The 
exterior derivative of (JT5j) is written by 

d/3 = d((?dp i )-d(F i dq i ). 

Then 



, a dG l , , dG\ , dF l , , <9F* , 
dp = — — ao 7 - A dpi + — — dpj A dpi — - — ag, A aoj — - — dp. A aoj. 
oqj dpj oqj dpj 

Here A is the wedge product. Equation (fTH|) can be presented in an equivalent form 



'dG j 9F\ , , 1 fdG j d&\ . , 1 /aP 
% <% / 2 V dp; <% / 2 V <9flj dqi 



Here we use the skew-symmetry of dqi A cig.,- and dpi A d!pj with respect index i and j . It is obvi- 
ous that conditions (fTTj) lead to the equation d(3 — 0. Equations (|T7jl are called the Helmholtz 
conditions jSHl EE HUE 122] for the phase space. 

Proposition 4. JTie dynamical system \1$ on the phase space R 2n , is a Hamiltonian sys- 
tem that is defined by the Hamiltonian H = H(q,p) if the form ilfy) is an exact form (3 = dH , 
where H = H(q,p) is a continuous differentiable unique function on the phase space. 

Proof. Suppose that the form (fT5j) is 

a JTT dH J dH J 
P = dH = ——dpi + -7— c%. 

dpi dqi 

Then the vector fields (G\ F l ) are 

dH dH 

&(q,p) = ir , F i (q,p) = -—. (19) 
dpi dq { 

If H = H(q,p) is a continuous differentiable function, then the conditions ()17|) are satisfied, 
and ([H]) is a Hamiltonian system. Substitution of (|T9"j) into (fT4"j) gives 

dqi_ _dH_ dpi _ _dH_ 

dt dpi ' dt dqi ' 

As the result, the equations of motion are uniquely defined by the Hamiltonian H. 



3.2 Fractional Hamiltonian Systems 

Fractional generalization of Hamiltonian systems has been suggested in Fractional ana- 
logue of the form ()15j) can be defined by 

A, = G\dpi) a - F\d qi ) a . (21) 
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Let us consider the equations of motion 

A t q t = G\q,p), B t p t = F l (q,p), (22) 

where A t and B t are the linear (or nonlinear) operators that have derivatives (integer or frac- 
tional order) with respect to time. As the simple examples, we can consider the total time 
derivatives A t = B t = d/dt and equation (|14j) . For the fractional derivatives A t = B t = D°, 

D? qi = G l (q,p), D? Pi = F\q,p). (23) 



Definition 4. The dynamical system \22]) on the phase space R 2n is called a fractional 
Hamiltonian system if Wfi) is a closed fractional form 

d a l3 a = 0, (24) 

where d a is the fractional exterior derivative. The system is called a fractional non-Hamiltonian 
system if \21)) is non-closed fractional form, i.e., d a (3 a ^ 0. 

The fractional exterior derivative for the phase space R 2n is defined as 

d a = {dq i ) a T>l + (dp i ) a T>l, a>0. (25) 
Let us consider a fractional generalization of the Helmholtz conditions. 

Proposition 5. If the right-hand sides of equations \22]) satisfy the conditions 

D" G i - D" G j = 0, D" G j + F l — 0, D" F i - D" F j = 0, (26) 

Pj Pi 'Hi Pj H3 Hi x ' 

then dynamical system \2§) is a fractional Hamiltonian system. 
Proof. This proposition has been proved in Ref. 



Proposition 6. The dynamical system \2fy) on the phase space R 2n is a fractional Hamil- 
tonian system with the Hamiltonian H = H(q,p) if \21\l is an exact fractional form 

(3 a = d a H, (27) 

where H = H(q,p) is a continuous differentiable function on the phase space. 

Proof. Suppose that the fractional form ()21|) is 

/3 a = d a H = (dp t ) a -D^H + (dq t ) a -DlH. 

Then 

(?(q,p) = T>°H, F i (q,p) = -T>ZH, 

and equation (jUJ) gives 

A t q i = B«H > B tPt = -r>lH. (28) 
These equations describe the motion of fractional Hamiltonian systems. 
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4 Hamilton's Approach 



4.1 Hamilton's Equations of Integer Order 

Let us consider Hamiltonian systems in the extended phase space M 2n+1 = R 1 x R n x R n 
of coordinates (t,q,p). The motion of systems is defined by stationary states of the action 
functional 

S[q,p] = J\pq-H(t,q,p)]dt, (29) 

where if is a Hamiltonian of the system, q = dq/dt, and both q and p are assumed to be 
independent functions of time. In classical mechanics, the trajectory of an object is derived by 
finding the path for which the action integral (J29)) is stationary (a minimum or a saddle point). 
In Hamilton's approach the action functional ()29|) can be written as 



S[q,p} = Juj hj (30) 

where 

u>h = pdq — Hdt. (31) 

The form (J3TJ) is called the Poincare-Cartan 1-form or the action 1-form. 

The Poincare-Cartan form looks like the integrand of the action or the Lagrangian. However, 
it is a differential form on the extended phase space M 2n+1 of (t, q,p), not a function. Once we 
integrate it over a curve C in M 2n+1 , we get the action: 

S[q,p] = [ uj h = [ B [pdq - H(t,q,p)dt}. (32) 

The integration is taken from A to B in the extended phase space M 2n+1 . 

Now suppose we integrate from A to B along two slightly different paths and take the 
difference to get a close loop integral. To evaluate this integral we can use Stokes' theorem 
In the language of differential forms, Stokes' theorem is written as 



to — I dui. (33) 

dM JM 

Here, M is a n-dimensional compact orientable manifold with boundary dM and u is a (n — In- 
form; 'd! is its exterior derivative. Note that M can be a submanifold of a larger space, so 
that Stokes' theorem actually implies a whole set of relations including the familiar Gauss and 
Stokes laws of ordinary vector calculus. 

Applying equation (J3*3*j) to the difference of actions computed along two neighboring paths 
with (q, t) fixed at the endpoints, we get 

5S[q,p}= [ du) h = [ (dp Adq- dH A dt) , (34) 

where a denotes the surface area in the extended phase space bounded by the two paths from 
A to B. 

The principle of stationary action states that 5S = for small variations about the true 
path, with (q, t) fixed at the end points. This will be true for arbitrary small variations, if and 
only if dujh = for the tangent vector along the extremal path. 
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We can consider the exterior derivative of Poincare-Cartan 1-form, and derive the equations 
of motion from the condition duh = 0. Using this condition, we get the Hamilton's equations 
of motion. This condition is equivalent to the stationary action principle 5S[q,p] = 0. 

Proposition 7. 

The exterior derivative of the Poincare-Cartan 1-form < f3Tj) is defined by the equation 

duo h = [D t p + D q H] dtAdq- [dq - D p Hdt] A dp, (35) 

where 

D ' = w °«=|< D -=w (36) 

Proof. The exterior derivative of the form (j3*Tj) can be calculated from the equation 

du h = d{pdq) — d(Hdt) = D t pdt Adq + D q pdq A dq + D p pdp A dq— 

- D t Hdt A dt - D q Hdq A dt - D p Hdp A dt. (37) 
Using dt A dt = 0, dp A dt = —dt A dp, and D q p = 0, we get 

duj h = [D t p + D q H] dtAdq - [D p pdq - D p Hdt] A dp. (38) 

The relation D p p = 1 gives 

du h = [D t p + D q H] dtAdq - [dq - D p Hdt] A dp. (39) 

Stationary Action Principle in Hamilton's Approach 

The trajectory of a Hamiltonian system can be derived by finding the path for which the Poincare- 
Cartan 1-form Uh is closed, i.e., 

du h = 0. (40) 
Using the stationary action principle (|4T)|) . we get the equations of motion 

dq - D p Hdt = 0, D t p = —D q H. (41) 

As the result, we obtain 



dq _ dp 
dt ~ p ' dt 
which are the well-known Hamilton's equations. 



^ = D P H, J = -D q H, (42) 



4.2 Fractional Hamilton's Equations 

The fractional generalization of the form (|31j) can be defined by 

u h , a = p(dq) a - H(dt) a . (43) 

Note that is a fractional 1-form that can be called a fractional Poincare-Cartan 1-form or 
fractional action 1-form. 



S 



We can consider the fractional exterior derivative of the form (JEJ), and use d a Uh, a = to 
obtain the fractional equations of motion. 



Proposition 8. The fractional exterior derivative of the fractional form (^3j) i 



is 



d a LU h>a = D?p + D"H (dt) a A(dq) a 



P 



l-a 



■(dq) a - D"H{dty 



A {dp) a . (44) 



_r(2-a) 

Proof. The fractional exterior derivative of the form (J4H|) is calculated by using the rule 

d k g 



k=0 



dx k 



and the relation 

For example, we have 
d a 



dx k 



[{dx) a ] =0 (k > 1). 



CO Qfc 

A\d Xl r\ = j2(d Xj r a Q p^w^ 

k=o ox j 



where 



(d Xj r a (dxr o (d« a') = (p^A*) (dx 3 r a (d Xl y 



KkJ r(l - a)T(k + 1) ' 

As the result, 

d a 0J hjOl = d a (p(dq) a ) - d a (H(dt) a ) = 

= (D?p)(dt) a A (dq) a + (D a q p)(dq) a A (dq) a + (D a p p)(dp) a A (dq) a - 

- (D?H)(dt) a A (dt) a - (D«H)(dq) a A (dt) a - (D a p H)(dp) a A (dt) a 

Here D a are the Riesz or Caputo fractional derivatives EHJ El E2]- Note the Riemann 
Liouville fractional derivative leads us to dependence of independent coordinates 



(45) 



D«p = pD«l^0. 



(46) 



Therefore the fractional equations are more complicated for Riemann- Liouville derivatives. 
Using (dt) a A (dt) a = 0, (dp) a A (dt) a = -(dt) a A (dp) a and D%p = for Riesz and Caputo 
derivatives, we can rewrite equation (J45)) in the form 

d a u Ka = \D?p + D°h] (dt) a A (dq) a - \(Dy)(dq) a - (D°H)(dt) a ] A (dp) a . (47) 



Substitution of 



into equation (j4Tj) gives (J33J). 



P 



l-a 



r(2-a) 



(46 
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Fractional Action Principle in Hamilton's Approach 

The trajectory of a dynamical system can be derived by finding the path for which the fractional 
Poincare-Cartan 1-form Uh, a is fractional closed form, i.e., 

d a u Ka = 0. (49) 

Here, we consider only fractional Hamiltonian systems. The non-Hamiltonian systems are 
considered in section 6. 

Using (jUJ) and flU, we get 

T{2- a) {dqT ~ D P H( < dt ^ = °' D ?P = -D%H. (50) 
As the result, we obtain 

(dq\ 
( dt 



T(2 - a)p a - l D«H, D^p = -D«H. (51) 



These equations are the fractional generalization of Hamilton's equations. 
For the fractional Poincare-Cartan 1-form 

uJ h , a =pP(dq) a -H(dt) a } (52) 

equations (Jol"j) are 

(dgV _ r(/3+l-q) ~ p _ 

[dt) ~ T(P + 1) P DpH > DtP -~ D « K (53) 

Note that we cannot use the rule of differentiating a composite functions for fractional deriva- 
tive DfpP . Therefore equations with (3^1 are more complicated than equation (|51|) . 



5 Lagrange's Approach 

5.1 Lagrange's Equations of Integer Order 

Suppose L(t,q,v) is a Lagrangian of dynamical system, where q is coordinate, and v is the 
velocity. Let us consider the variational problem on the extremum of the action functional 

S [q,v] = J L(t,q,v)dt, (54) 

under the additional condition 

q = v, (55) 

where both q and v are assumed to be independent functions of time. In this case, p play 
the role of independent Lagrange multipliers. Obviously, the indicated variational problem is 
equivalent to the problem on the extremum of the action, 

S[q, v,p] = J [L(t, q, v) + p(q - v))dt, (56) 
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where already all the variables q, v, p have to be varied. The corresponding Lagrange's equations 
are 

q = v, V = Tv (57) 

We can introduce the extended Hamiltonian in the space of variables (t,q,p,v) as 

H*(t,q,p,v) =pv - L(t,q,v). (58) 
The corresponding extended Poincare-Cartan 1-form is defined by 

uJh* = pdq — H*dt = pdq + Ldt — pvdt. (59) 

Proposition 9. 

The exterior derivative of the form is defined by the equation 

duJh* = [D t p — D q L}dt Adq — [dq — vdt] A dp — [p — D v L]dv A dt. (60) 

Proof. The exterior derivative of (|59|) is 

dcdh* — d(pdq) — d(H*dt) = D t pdt A dq + D q pdq A dq + D p pdp A dq + D v pdv A dq— 

- D t H*dt Adt- D q H*dq A dt - D p H*dp A dt - D v H*dv A dt. (61) 
Using dt Adt = 0, dq Adt = —dt A dq, and D q p = 0, D v p = 0, equation (joT|) gives 

duj h * = [D t p + D q H*}dt Adq- [D p pdq - D p H*dt] A dp - D v H*dv A dt. (62) 

From (j5HJ), we get 

D q H* = D q [pv — L) — —D q L, 
D P H* = D p [pv — L] = vD p p — D p L(t, q, v) = vD p p = v, 
D V H* = D v [pv — L] — pD v v — D V L = p — D V L. 
As the result, we obtain (loT7|) . 

Stationary Action Principle in Lagrange's Approach 

The trajectory of a dynamical system can be derived by finding the path for which the form A5y\) 
is closed , i.e., 

duJh* = 0. (63) 

From equations (JoPj) and (JBMj) . we get 

D t p - D q L = 0, dq- vdt = 0, p - D V L = 0. (64) 

It is easy to see that equation (|64|) coincides with the Lagrange's equations (|57|) that can be 
presented as 

' d 

-D V L 



DqL — 



dt 



= 0. (65) 

v=q 



As the result, we obtain 

dL d ( dL\ , . 

which are the Euler-Lagrange equations for Lagrangian system. 
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5.2 Fractional Lagrange's Equations 

Suppose L(t, q, v) is a Lagrangian of the system, and the extended Hamiltonian is defined by 



H*(t,q,p,v) = pv p - L(t,q,v), 
where (3 is a positive power. Let us define 

u h *a = p(dq) a - H*(dt) a = p(dq) a + L{dt) a - pv?{dt) a , 
which is a fractional generalization of the extended Poincare-Cartan 1-form ()59|). 

Proposition 10. 

The fractional exterior derivative of the fractional 1-form Ab'ty) is defined by 



(67) 



(68) 



d a uJh* 



D°p — D°L (dt) a A(dq) a - 



- (Dy) \(dq) a - V P(dt) a ] A (dp) a - \pDyP - D«L] (dv) a A (dt) c 



(69) 



Proof. The fractional exterior derivative of the form (|68|) gives 

d a u h * a = d a [pd a q] - d[H*(dt) a ] = (D?p)(dt) a A (dq) a + (D°p)(dq) a A (dq) a + 

+ {D a p p){dp) a A (dq) a + {D a v p){dv) a A (dq) a - 
-(D?H*)(dt) a A {dt) a - (D a q H*)(dq) a A (dt) a - 

- (D^H*)(dp) a A (dt) a - (D*H*)(dv) a A (dt) a . (70) 

Using (dt) a A {dt) a = 0, {dq) a A {dt) a = -{dt) a A {dq) a , and D q p = 0, D v p = for Riesz or 
Caputo fractional derivatives, we can rewrite (|68|) as 

d a w h * a = [Dfp + D«H*] {dt) a A {dq) a - 

- \D«p){dq) a - [D*H*){dt) a ] A {dp) a - {D«H*){dv) a A {dt) a . (71) 
Using the extended Hamiltonian (JBTj) and properties of Riesz and Caputo derivatives, we have 

D«H* = D«\pvP-L] = -D«L, 

DpH* = D«\pv? — L] = v?D a p p - DpL(t, q, v) = v?D«p, 
D^H* = D^lpv? — L] = pD%v p - D a u L. 
As the result, we obtain (JBTJj) . 

Fractional Action Principle in Lagrange's Approach 

The trajectory of fractional dynamical systems can be derived by finding the path for which the 
fractional extended Poincare-Cartan 1-form \68\) is closed form, i.e., 



d a uih* 



0. 



(72) 
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Using (jfinj) and |J75|). we have 

p t a p) - -D^L = 0, - ?/(dt) a = 0, V D a y - D a v L = 0. 

From the relation 



„ r(,3 + i) 

" r(/3+l-a) ' 



where (3 > — 1, we get 



Substitution of the third equation from (|75j) into the first one gives 



rpg + i-q) 
r(/? + i) A 



v = (q) a ^' 



Equation (fTUj) has the dependence 

vf = (q) a . 

It is easy to see that equation (|7H|) looks unusually even for /3 = 1. Therefore we use /3 
the Hamiltonian f)67j) and the form (|68p. 
Using /3 = a in equations flHZJ), (JSHJ), and 



r(2-a)' 

we have the fractional extended Lagrange's equations 

(D?p) = D%L, ^ = ||p p = T(2-a)D«L. 

Substituting the third equation from (jHUj) into the first one, we obtain 

D«L - r(2 - a)D? [L>«L]„ = . = 0, 

that is fractional Euler- Lagrange equations. As the result, the fractional equations of 
for Lagrangian systems are presented by 

E?L(t,q,q) = 0, z=l,...,n, 

where 

E? = Dl-T{2-a)D«Dl. 
For a = 1, equations (|%T|) are the Euler-Lagrange equations (Jfifij) . 
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6 Fractional Non-Hamiltonian Systems 



6.1 Helmholtz Conditions 

In this subsection, we consider the brief review of Helmholtz conditions [36, 37 to fix notation 
and provide a convenient reference. The Helmholtz equations give the necessary and sufficient 
conditions for equations to be the Euler-Lagrange equations that can be derived from station- 
ary action principle. 

Proposition 11. The necessary and sufficient conditions for 

F i (t,q,q,...,qW) = o i i = l,..., n (83) 
to be equations that can be derived from the stationary action principle are 

where qf' = d k qi/dt k , i,j = 1, ...n and 

(k x = k\ d B " (fc) d k 

[m) m\(k-m)V dt dt dq\ k - 1Y 
Proof. This proposition is proved in |37j . 

For simple example, let us consider the equations 

F t (t,q,q) = 0, (87) 

Conditions ()84jl and (|85jl have the form 

dqj dqi dt \Btji J 
BF- BF 

^ + ^ = 0, i,j = l,...,n, (89) 



ddj 8q { 



where 



Equations (|8"9j) gives 



d B . 8 , nnS 

Tt = 8i + qi ^ (90) 



B 2 Fi 

0, i,j,K = l,...,n. (91) 



BqjBq K 

These conditions are satisfied for the linear dependence Fj with respect to q, i.e., 

Fi = Cij(t, q)<jj + Di(t, q) = 0, z = l,...,ra. (92) 
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Corollary 1. The necessary and sufficient conditions to derive equations UJ^) from the 
stationary action principle have the form 

Cij = —Cji, (93) 

dCjrj dC, K dC K , , , 

dq K dq t dqj 

at oqj oqi 

Proof. Substitution of equation into equations (jHHI) and (jHH) yelds flg3j) , (l9"3j) . and (1931). 

Corollary 2. The necessary and sufficient conditions for 

Fi(t,q,q,q) = 0, i = l,...,n (96) 
to fre equations that can be derived from stationary action principle are 

dFi dF,; . . 

dqj + % 2 dt (sg'j) °' ^ 

dF\_dF) d_ (9Fj\ _ <f_ fdFj\ () fg9 N 
9gj (it \ 9gj y (it 2 y <9g'j y 

where 

d d . d .. d 
dt dt ^ % dq, 1 dqi 
Note that using equation (|9*Hj) condition ()99j) can be rewritten in the more symmetric form 



dFi dF, 1 d (dF, dF, 



dqj dqi 2 dt \ dqj dq\ 



0. (101) 



6.2 Non-Lagrangian Systems 



Definition 5. A dynamical system is called non-Lagrangian system if the equations of motion 
\8'J\) cannot be represented in the form 

a-^p* « m)=0 - (102) 

wt/i some function L = L(t, q, q, q^) , where = d k q/dt k . 

It is well-known that the equations of second order cannot be presented as 

EiL(t,q,q) = 0, i = l,...,n, (103) 
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where E{ is the Euler-Lagrange operator 



E^D qi -D t D qi . (104) 

In the general case, the Lagrange's equations have the additional term Qi(t,q,q) which is a 
generalized non-potential force. This force cannot be presented as Qi = EJJ for some function 
U = U(t,q,q). In general, the Euler-Lagrange equations |H] are 

E i L(t,q,q)+Q i = 0. (105) 

If we consider non-potential forces and non-Lagrangian systems, then the nonholonomic 
variational equation suggested by L.I. Sedov jlOJ HU H21 EE] should be used instead of stationary 
action principle. 

6.3 Non-Hamiltonian Systems and Friction Force 

In general, the phase-space equations of motion cannot be presented in the form 

ft = D Pi H, Pi = —D qi H, (106) 

where H = H(t, q,p) is a smooth function. The Hamilton's equations are written as 

qi = D Pi H + G%q,p), p i = -D 9i H + F i (t,q,p), (107) 

where H = H(t,q,p) is a Hamiltonian of the system. For example, H(t,q,p) = T(p) + U(q), 
where T(p) is a kinetic energy, and U(q) is potential energy of the system. The functions 
G l (t,q,p) and F % {t, q,p) describe the non-potential forces which act on the system. For me- 
chanical systems, we can consider G l (t,q,p) = 0. If the functions G l (t,q,p) and F l (t,q,p) do 
not satisfy the Helmholtz conditions 1)17)1. then ()107)) is a non-Hamiltonian system. 

In general, the exterior derivative of the Poincare-Cartan 1-form is not equal to zero (duh ^ 
0). This derivative is equal to differential 2-form 9 that is defined by non-potential forces 

6 = F i (t,q,p)dtAdq i -G i (t,q,p)dtAdp i (108) 

for the non-Hamiltonian system (J107j) . For example, the linear friction force F % = —^ipi gives 

9 = -j Pi dt A dq { . (109) 

Proposition 12. The differential 2-form 9 of non-potential forces is non-closed form. 

Proof. If differential 2-form 9 is a closed form (d9 = 0) on a contractible open subset X of 
R 2n , then the form is the exact form such that a function h = h(t, q,p) exists, and 9 = dh. In 
this case, we have a new Poincare-Cartan 1-form 

u' h = uj h + h, 

such that duj' = 0, and the system is Hamiltonian. 
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There is a generalization of stationary action principle for the systems with non-potential 
forces. 

Action Principle for non-Hamiltonian Systems 

The trajectory of a non-Hamiltonian system can be derived by finding the path for which the 
exterior derivative of the action 1-form ( fffl)) is equal to the non-closed 2-form MQty) . i.e., 

duj h = 6. (110) 

Equations (|35|h (jl(J8|) and (jll(J|) give the equations of motion (ll(J7|) for non-Hamiltonian 

system. 

6.4 Fractional Generalization of non-Hamiltonian Systems 

Let us define a fractional generalization of the form (|1(J8|) by 

9 a = F l (t,q,p)(dt) a A (d gi ) a - G\t,q,p){dtT A (d Pi ) a . (Ill) 
This form allows us to derive fractional equations of motion for non-Hamiltonian systems. 

Fractional Action Principle for non-Hamiltonian Systems 

The trajectory of a fractional system subjected by non-potential forces can be derived by finding 
the path for which the fractional exterior derivative of the fractional action 1-form ^Bjj is equal 
to non- closed fractional 2-form ill!]) , i.e. 

d a u ha = 6 a . (112) 

Using (jSD, (HTTP and (fTT2D . we get 

T ^* a) {dq) a -D a p H{dty = G(t,q,p)(dt) a , D?p = -D a q H + F(t,q,p). (113) 
As the result, we obtain 



dqi 
dt 



r(2 - a)pr X D a v M + G% q,p), D? Pi = -D«H + F% q,p). (114) 



These equations are the fractional generalization of equations of motion for non-Hamiltonian 
systems. 
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7 Conclusion 



In this paper, we define a fractional exterior derivative for calculus of variations. Hamilto- 
nian and Lagrangian approaches are considered. Hamilton's and Lagrange's equations with 
fractional derivatives are derived from the stationary action principles. We prove that frac- 
tional equations can be derived from action which has only integer derivatives. Derivatives of 
noninteger order appear by the fractional variation of Lagrangian and Hamiltonian. 

Application of fractional variational calculus can be connected with a generalization of 
variational problems. The gradient systems (GS) form a restricted class of ordinary differential 
equations. Equations for GS can be defined by one function - potential. Therefore the study of 
GS can be reduced to research of potential. As a physical example, the ways of some chemical 
reactions are defined from the analysis of potential energy surfaces jlSJ The fractional 
gradient systems (FGS) have been suggested in [231- It was proved that GS are a special 
case of such systems. FGS includes a wide class of non-gradient systems. For example, the 
Lorenz and Rossler equations are fractional gradient systems (231- Therefore the study of the 
non-gradient system which are FGS can be reduced to research of potential. 

We can assume that the ways of some chemical reactions with dissipation and systems with 
deterministic chaos can be considered by the analysis of fractional potential surfaces. Let us note 
the interesting property of potential surfaces for systems with strange attractors. The surfaces 
of the stationary states of the Lorenz and Rossler equations separate the three-dimensional 
Euclidean space into some number of areas (23] • We have eight areas for the Lorenz equations 
and four areas for the Rossler equations. This separation has the interesting property: all 
regions are connected with each other [23]. Beginning movement from one of the areas, it is 
possible to appear in any other area, not crossing a surface. Any two points from different 
areas can be connected by a curve which does not cross a surface. 

The fractional variations can be used to define the fractional generalization of gradient type 
equations that have the wide application for the description dissipative structures |I7| 138] . 
The fractional gradient type equations are generalization of FGS [23] from ordinary differential 
equations into partial differential equations. We plan to realize this generalization in the next 
paper by using de Donder-Weyl Hamiltonian and Poincare-Cartan n-form. 
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